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We prove that the hierarchical fractal model recently proposed for describing the stretched 
polymers [A. N. Samukhin et al, Phys. Rev. Lett. 78, 326(1997)] is equivalent to a one-dimensional 
chain with hierarchical aperiodic structure. By use of the transfer matrix technique we calculate 
the electronic transmission and the dc conductance. We find that there exist sharp-edged trans- 
mission subbands and gaps, but the transmission subbands are rich in substructures which show 
self-similarity. The temperature dependence of conductance g{T) is sensitive to the variation of the 
Fermi level and to the structure parameters m and n of the original fractal structure. The relation- 
ship of the obtained results to the experimental data and other theoretical works is discussed. 
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I. INTRODUCTION 

It is widely accepted that polymers represent complicated structures which consist of coupled one-dimensional 
(ID) chains, and many unusual properties can be understood based on their quasi-lD feature. During last decades 
analytical and numerical studies revealed that the phenomenon of Anderson localization is a characteristic feature of 
ID systems. For a quasi-lD system the interchain coupling tends to delocalize the electronic states, which may lead 
to extended states when interchain diffusion rate cu exceeds a threshold value ujja- Moreover, conductance of several 
polymers, such as stretched polyacetyllene and polyanilline, is enhanced by a few orders of magnitude upon doping 
and can even reach high values, typical of metals, while the temperature dependence of static conductance (gdc) widely 
vary from sample to sampld3~a. 

In order to understand the effect of chain_windings and random interchain couplings on the electronic transport 
properties in quasi-lD polymers, Xiong et aBx3 have investigated the transmission coefficient and conductance for 
models of winding and randomly coupled polymers. It is proved that for such systems there exist conducting subbands 
with "dilute" structure which can account for the high values jof conductance and its insulator-like temperature 
dependence of some polymers. On the other hand, Samukhin et. aOsng gested that the peculiar features of conducting 
polymers can be derived based on a special hierarchical structure-a nearly ID fractal. ■— . 

In this paper we investigate the transmission spectrum of the hierarchical structure proposed irJlil. We prove 
that this fractal network is equivalent to a ID chain with hierarchical aperiodic structure for the transport of electrons 
in the stretched direction. By solving the Schrodinger equation of such an effective ID model by the standard transfer 
matrix technique, we obtain the transmission coefficient as a function of energy for diffscent stages (or size) of the 
structure. We find that there also exist "dilute" conducting subbands as those in RefB andllJ, but the structure 
in these subbands exhibits some kind of self-similarity and become more complicated when the stage of the system 
increases. As a result, the calculated conductance sensitively depends on the structure parameters (m, n) of the 
original hierarchical fractal and on the position of the Fermi level. 

The paper is organized as follows: In the next section we present the derivation of the equivalent ID model and 
the basic formalism. In Section III the calculated transmission spectrum and conductance are illustrated. Section IV 
is devoted to a brief summary and discussion. 

II. BASIC FORMALISM 

According to Samukhin et. al, the hierarchical model for the stretched polymer is constructed in the following 
way: as a first stage, n primitive bonds are taken to form an ri-length chain along stretched direction, and m threads of 
such chains are combined at both ends to form an (m, n) boundle with two common terminals; in the following stages 
the boundlcs formed in the previous stage are taken as the elementary bonds and the same procedure is performed 
to form new boundlcs. We are interested in the variation of the transmission spectrum in increasing the number of 
stages N, i.e. the size of the system. In measuring the conductance one needs to link electrodes to both ends, (see 
Fig.l). We use two integers, to label the sites, where i indicates the position in the longitudinal direction and 
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j in the transverse direction. If the number of stages of the system is TV, the longitudinal index i is varying in the 
range of [0,71^]. For a given i, the number of sites in the transverse direction is 

, if mod(i,n) 7^ 0; 

PN{i) = { rn^-""(*) if mod(i, = and mod(i, n""(')+i) ^ 0,aN{i) is integer in (0, N); (1) 

1, if i < or i > , 

where mod(a;, y) denotes the remainder on the division of x by y. Thus, in the site index (i, j), j is ranging in [1, pff[i)]. 
From the tight binding approach, the Schrodinger equation can be written as 

Ei>i,j = iijipij + ^ itj,&{t,j)i^s(i,j)^ (2) 

where V'ij is the coefficient of the wave function on site S{i,j) denotes the position of next neighbors of (i, j), 

eij is its site energy, the hopping between sites and d{i,j). As we are only interested in the topological 

structure of the system, we can assume that the site energies and the nearest-neighbor hoppings are independent of 
the position. In this case one can set the site energy to be the energy origin and the hopping to be the energy unit. 
By summing over the transverse sites for a given i, Eq. (|^) becomes 

E'ipi = /3Ar(i,i + l)'ipi+i + (iNihi- (3) 

where 

, PN{i) 



VPnW ,=i 



and 

I PNii) PN{i') 



From Eq. p) it is easy to see that 



PN{i') ' V Pn{^) 



/3Ar(«,i') = (4) 

where aAr(j) is the largest power of factor n""'*-' of integer i in the range [1, — 1], as defined in Eq. (|l|), and we 
set a^ii) = iV for i < and i > . 

FromEq. (|) one can see that the fractal structure reduces to a ID effective Hamiltonian with special aperiodic 
structure. In the derivation we have adopted the summation over the transverse sites which erases the degree of 
freedom of electrons in this direction. This has no effect on the longitudinal transport properties for the one-channel 
incoming and outgoing leads of the present geometry. The situation becomes different if multichannels are connected 
to the system or there exists on-site or off-diagonal disorder in the system which can produce an energy-dependent 
terms in the reduced Hamiltonian. In this paper we only focus on the present geometry and the multichannel and 
disordered cases will be considered in separated ones. 

The transmission coefficient in such a one-channel structure can be calculated by the standard transfer matrix 
technique. The transfer matrix can be defined as follows 

, \ / ^ °w(i)-aw(' + l)l l°jv(')-°iv('-l)|-|°jv(')-°w('+l)l \ / ; 

Vi+1 \ _ / Em 2 —m 2 \ / iPi 



V'J / V 1 J 



So along the ID chain with stage number N, the wave function of an electron in two leads are related to each other 
via a total transfer matrix 



z=0 
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We demonstrate such an effective ID chain structure in Fig. 2. This is a one-to-one mapping from the original 
network of Fig.l, and the resuhant ID chain shares the hierarchical structure at its longitudinal direction and the 
mirror reflecting symmetry as well. The self-similar features of the original structure are also kept in the mapping. 

The pAr(i)-structure in the longitudinal direction of the original system can be described by the following series: 



1 iv iv — i I\ Jy Jy I\ — L I\ 

l,m ,---,m ,m ,m ,■ ■ ■ ,m ,■ ■ ■ ,m ,m ,---,m ,m (7) 
^ , " , " , " 

n— 1 n—1 n—1 

~- V ' 

n-1 

One can see that this is a embedded multi-period structure. As a consequence the hoppings in the effective ID chain 
are arranged as in the following series: 

1, 1, • • • , 1, , 1, • • • , 1, m^, • • • , m^, 1, • • • , 1, , 1, • • • , 1, mi , ■■ ■ . (8) 

n—2 n—2 n—2 n—2 



n-2 

Thus, one has the following recursion relation of the transfer matrix for n > 2: 

T(iV) = BM{n^)BN{n^ - l)BN-i{n^-^ - 1)C{N - l)[C{N - 1) 

X F{N - 1)]"-2C(7V - l)B^i_,(l)i?jv(l)Bjv(0), (9) 

where 

C{N - 1) = i?^i_,(n^-i)T(iV - l)B],l,{0), 

and 



F{N -1) = V^^^ ^ ^ 



One can find that the effective ID chain contains bonds among which there are (n — 2)n^ ^ bonds of 
hopping 1, 2 bonds of hopping V , and 2n^~^~^(n — 1) bonds of hoppings Vm* for i = 1, 2, • • • , — 1. These 
bonds are arranged in a special aperiodic way so that one may expect a self-similar and Cantor-set-like transmission 
spectrum. Prom the transfer matrix T, we can easily calculate the transmission coefiicient for electrons from the left 
to the right lead: 

4sin^A; 

' " [T21 - T12 + (T22 - Tn) cos fc]2 + (Tn + Taa)^ sin^ k ' ^ ' 

where k = arccos(£'/2) is the wave vector of wave function in the leads. Thus, the conductance of the original network 
at the longitudinal direction is obtained according to Landauer formula 

5(£;) = ^ fr{E)^dE (11) 



where / is the Fermi distribution function 



III. TRANSMISSION SPECTRUM AND TEMPERATURE DEPENDENCE OF CONDUCTANCE 

By using the above formulas we calculate the transmission coefficient as a function of energy for different 
structure parameters (m, n) and different stage number. Typical results of transmission are illustrated in Fig. 3 and 
Fig. 4. We can see that when the stage number of the structure is low (say < 2), t{E) exhibits quite continuous 
variation with broad peaks and dips. However, when the stage number A > 3 the whole energy interval [-2,2] is 
divided into "subbands" which allow electron transmission and "gaps" of nearly zero transmission with sharp edges. 
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By increasing N, the subbands are further divided into "minibands" and "minigaps", but the edges of the main 
subbands are almost independent of TV. For large N, there are a large number of sparsely distributed peaks with t(E) 
equal to or near unity which indicates almost complete transition corresponding to unscattered electronic states. This 
structure of the transmission spectrum also shows the self-similarity and Cantor-set-like behavior, like the situation 
in the Fibbonaci series and other quasiperiodic or aperiodic systems. By the comparison between Fig. 3 and Fig. 4, 
we find that although the number of the main subbands is the same but the width of the central main subband is 
slightly widened and the other two are shifted away from the band center by increasing n, which changes the number 
of the repeated bonds in the series. In the inset of Fig. 1(a) we also show the TV dependence of the transmission 
coefficient for energy E — 1.57 which lies in the main band. It can be seen that for this energy by increasing N the 
transmission first decreases but then increases. This behavior is sensitive to the position of E as can be seen from the 
complicated structure in the transmission spectrum for large N. 

For the polymer networks, an interesting problem is to investigate the scaling behavior and the temperature 
dependence of the dc conductance. Since the transmission subbands exhibit "dilute" characteristics with a lot of 
transmission holes within them, the temperature dependence of the conductance is sensitive to the position of the 
Fermi level. From Eq. (|l^), for low temperature the conductance as a function of the Fermi level shows almost the 
same characteristics as those for the transmission coefficient as a function of energy. From Figs. 3 and 4 one can see 
that there exist energy points for which the transmission of large system is greater than that of small system. As a 
consequence for the system with Fermi level at these points the conductance increases with the system size provided 
that the temperature is low enough. This may be regarded as the "metallic" scaling behavior. There are still many 
other energy points for which the scaling behavior of the conductance is "insulating" as it decreases with the system 
size. However, the energy regimes for "metallic" and "insulating" states are strongly mixed with each other and form 
a complicated Cantor-set-like structure. If the temperature becomes higher, such a structure is smoothened by the 
thermal averaging. In Figs. 5 and 6 we plot the Fermi level dependence of the conductance at temperature T = 400K 
for different structures and different system size. It can be seen that at this temperature the curves become smooth 
and in the whole energy range the conductance decreases with the increase of the system size, showing the insulating 
behavior. This is because more transmission holes appear in the transmission subbands when N increases and this 
reduces the average of the transmission. In Fig. 7 we plot the temperature dependence of conductance g(T) for one 
structure but with different Fermi levels. The curves manifest that the conductance can increase or decrease with 
the decrease of the temperature, strongly depending on the Fermi level Ef. The orders of magnitude of the variation 
of conductance at the range of [50K, 400K] could vary from 1 to 4. Only for carefully chosen values of Ef (usually 
near edges of subbands for a fflW^n structure (m,n)), one could get a curve of logg(r) - ^hose slope can be 

fitted to the experimental dataoa. In fact, the realistic stretched polymers should contain tremendous number of such 
networks in different stages, and doping might adjust Fermi level as well as the distribution of (m, n). This may cause 
the deviation of the conductance calculated results in this simple model from the observed data. Nevertheless, as one 
can see from Fig. 7, our calculation displays a curve of good fitting to expression Ing^c oc T~^/^ within quite large 
range of temperature, while the solution of the model from the critical percolation approach in Ref.B demonstrates 

that gdc oc °r^i^ — which slightly diverges from experimental results. 



IV. CONCLUSIONS AND DISCUSSION 



In summary, we have shown that the hierarchical fractal network of Fig. 1 is equivalent to a ID aperiodic 
series if only one-channel transport properties are under consideration. This equivalence is derived by summing over 
the transverse degree of freedom and, as a consequence, does not imply the one-dimensionality of this structure. The 
self-similarity and the aperiodicity of the original network are kept in this mapping. We calculate the transmission 
spectrum and temperature dependence of the conductance for systems with different structure parameters and various 
size. The transmission spectrum manifests Cantor-set-like characteristics, with complicated structures made of mixed 
transmission peaks and gaps, like most of the ID quasiperiodic and aperiodic systems. Such a structure is also 
reflected in the mixed insulating and metallic scaling behavior of the conductance at low temperature. At higher 
temperature such a structure in the conductance is smoothened by the thermal averaging and the scaling behavior is 
insulator like in the whole energy range. This implies some type of metal-insulator transition at special values of the 
Fermi level, in agreement with the analysis in Ref.E2l. The magnitude and temperature dependence of the conductance 
strongly depend on the position of the Fermi level. By adjusting the value of Ef, one can reproduce the temperature 
dependence observed in the experiments. This simplified model can catch several features of polymer qualitatively 
but a complete description needs more sophisticated models. 
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FIG. 1. Hierarchical fractal network with structure parameters m = 3, n = 4 and stage number (a) A'' = 1, (b) N — 2. 

FIG. 2. The equivalent ID chain for the hierarchical structure shown iu Fig. 1 with the same parameters. The numbers 
below the lines indicate the index of sites in the longitudinal direction, and the rmmbers above the lines arc the effective hopping 
integral of the corresponding bonds. 

FIG. 3. The transmission coefRcient as a function of energy for structures with m = 2, n = 5, and (a) A'^ = 1, (b) N = 2, 
(c) N ^ 3, {d) N = 4, (e) N = 5. The inset of (a) shows the N dependence of r(.B) for E = 1.57. 

FIG. 4. The transmission coefRcient as a function of energy for structures with m = 2, n = 6, and (a) N = 1, (h) N = 2, 
(c) N = 3, (d) N ^ 4. 

FIG. 5. The conductance as a function of Fermi level at temperature T = 400K for structures with m = 2, n = 5, and (a) 
N = 1, (h) N = 2, (c) N = 3, {d) N = 4, (e) N = 5. The hopping integral of the original network is set to be leV. 

FIG. 6. The conductance as a function of Fermi level at temperature T = 400K for structures with m = 2, n = 6, and (a) 
AT = 1, (b) AT = 2, (c) AT = 3, (d) N = 4. The hopping integral of the original network is set to be leV. 

FIG. 7. The conductance as a function of temperature of system with structure parameters m = 2, n = 5 and AT = 5 for 
different Fermi levels. 
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